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The probability distribution of the total transmission is studied for waves multiple scattered from 
a random, static configuration of scatterers. A theoretical study of the second and third cumulant of 
this distribution is presented. Within a diagrammatic approach a theory is developed which relates 
the third cumulant normalized to the average, {(T^)), to the normalized second cumulant {(T^)). 
For a broad Gaussian beam profile it is found that ((T^)) — ^{(Ta))'^. This is in good agreement 
with data of optical experiments. 
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Q ■ I. INTRODUCTION 

o 

CO ' Multiple scattering in disordered systems is a field of wide interest; it is studied in electronic, microwave and 
optical systems. In the multiple scattering regime the main transport in transmission is through diffusion. Yet 
. interference processes, possible by the underlying wave character, play an important role. This interference leads 

" to interesting effects such as the enhanced backscatter cone [|l|,||,D, short-, and long-range correlations [||j|J|,0,||, 

■ and strong localization. If interference occurs between diffusion paths, it causes large ffuctuations. Most famous 
] are the sample-to-sample fluctuations in the conductance of electronic systems, the so called universal conductance 

C\J ■ fluctuations but also other transmission quantities are influenced by interference. 

Recently, the attention is not only drawn to the variance of the fluctuations, but to the entire distribution functions. 
\l Examples are the intensity distribution in speckle patterns for classical waves | [l^ , [l3t and the conductance distribution 
for electronic systems [^dj. The size of the fluctuations and the shape of the distribution is related to the 'distance' 
from the localization transition. Far from localization, diffusion channels are almost uncorrelated and fluctuations 
are small (except for the optical speckle pattern in the angular resolved transmission). The correlation between the 
channels increases if the localization is approached. The relevant parameter is the inverse dimensionless conductance 

■ l/g, which can be interpreted as the chance that two channels interfere. The dimensionless conductance can be 
(~| ' expressed in the thickness of the sample L, the mean free path £, and the number of channels N, 
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The number of channels is calculated in analogy with a waveguide, where it is unambigiously deflned. In the diffuse 
^ ' mesoscopic regime will be a small parameter of our perturbation theory (experimentally this proved fully justifled, 
5^ , as there g ^ 10'^, p5|). Close to Anderson localization g approaches unity, and fluctuations increase. The central 
question is how the distributions change as the strong localization regime is approached . 

Let us briefly review some characteristics of optical transmission distributions functions in the regime of moderate 
g. In the study of mesoscopic systems using light scatterering, one takes a small sample with static scatterers. In 
order to average, one needs to sum over a large number of scatterer conflgurations, in practice this is done by varying 
an external parameter such that the interference pattern is completely changed. Whereas in electronic systems it 
is common to vary the Fermi energy or apply a magnetic field, in optics one usually varies the wavelength of the 
light. In contrast to electronic systems, not only the conductance but three different transmission quantities can be 
measured in optical systems. (An exception is the very recent observation of electron speckle by Gao et al. p7|.) 
First, the angular resolved transmission can be considered: If a laser shines on the sample, a speckle pattern is seen in 
the transmission. If one measures the intensity in an outgoing direction b, this corresponds to measuring the angular 
transmission coefficient Tab, where a denotes the incoming and b the outgoing channel. The intensities in the speckle 
pattern have in zero'th order in g an exponential distribution 

P{Tab) = -^exp(-^), (2) 



(Tab) V (Tab) 
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which is also known as the Rayleigh law. Deviations from this law occur if interference between the transmitted 
intensities, the diffusons, is taken into account: the higher moments of the distribution function increase. For large 
values of the transmission a stretched exponential was observed which was recently predicted as PiTab) oc 
exp(-2^gr,fa/(T,fa) ). 

One obtains another transmission quantity, if one collects all the outgoing light, and also all incoming directions 
are used, i.e. using a diffuse incoming waves instead of a plane wave. In this case one summes over all incoming and 
outgoing channels and the conductance of the sample, g = ^ Tab is measured. It corresponds to the well known 
conductance measurement for electron systems. The conductance has roughly a Gaussian distribution. The absolute 
variance is larger than expected classically because of interference and it does not depend on the sample parameters 
(hence the name universal conductance fluctuations) Equivalently, the relative variance is proportional to 

g~'^ ■ The full distribution function of the conductance was studied by Altshuler et al. [Q, who predict a log- normal 
tail for the tail of the distribution function. For the lower cumulants they predict that (5")cum c>c (5^^"). Note, 
however, that the prefactors maybe zero, as Macedo found that {g^)cum fx 

We are, however, interested in yet another transmission quantity. It is obtained by spatially integrating the speckle 
pattern on the outgoing side, and taking an incoming plane wave. One thus obtains an intermediate quantity between 
the angular transmission Tab and the conductance X^ab-^a^' -"-^ termed the total transmission Ta = '^^'^ab- This 
quantity is the subject of this paper. The total transmission is a constant superposed with fluctuations. In first order 
of g~^ the fluctuations have a Gaussian distribution [slJT^. The relative variance of this distribution is proportional 
to g~^, it is thus a factor g larger than for the conductance fluctuations. This sensitivity of the total transmission 
to interference processes and its simple limiting behavior (as compared to the angular resolved transmission) make it 
an ideal quantity to study mesoscopic transmission. Its full distribution was studied numerically in 2D by Edrei it et 
al. who recovered the Gaussian distribution function, tending towards a log-normal distribution near the Anderson 
transition [|l9|. The full distribution function was recently derived [|6|. It was shown that it has a log- normal 
distribution growth, and an exponential tail. 

Recently the third cumulant of the distribution was found experimentally by De Boer et al. [p^ . In this paper 
we present the theoretical details of that work. We focus on the Gaussian distribution and the deviation from the 
Gaussian due to the presence of the third cumulant. The structure of this paper is as follows. The diffusion in optical 
CT^tems is described in section ^ in section III the character of the probability distribution is discussed. In section 
IV and |v|+[vl, respectively, we calculate the second and third c umula nt of the probability distribution. Next we 
calculate experimental corrections to our result in section VU and VIII, after which we compare our results with the 
experimental data in section IX. We close with a discussion. 



II. DIFFUSE TRANSPORT 



Consider the transport of light through a three dimensional slab with static scatterers in a random configuration. 
The slab has a thickness L. Diffuse mesoscopic regime is characterized by A ^ ^ ^ L. Here £ denotes the mean free 
path and A denotes the wavelength in the medium. Isotropic scattering is assumed. We also work in the scalar wave 
approximation. It is known that the two independent polarization directions of light effectively double the number 
of channels N in the problem as compared to the scalar wave case . Apart from a small correction term, our final 
result does not depend on the number of channels and the doubling. 

We take the depth in the the slab as the z-coordinate, the slab thus corresponds to < 2 < L. A plane scalar wave 
with unit flux and area A impinges on the sample from direction a. It is given by 

V'r„(r) = rr-i eyipiiPaR + ik^iaz), z < 0, (3) 

^Ak^a 

where k is the wavenumber, r = (a;, y) is the transversal coordinate. Pa is the transverse component of the momentum, 
and /ia = ^1 — Pq/A^ = cos 6*0, where 9a is the angle with respect to the z— axis. The number of channels depends 
on the area of the incoming beam. Dividing the total area of the beam with diameter po in small coherence regions of 
area A^, one obtains in analogy with waveguides N = 2 x k^p^/A, or the Weyl formula N = 2 x k'^A/A-K. The factor 
2 comes from the two independent polarization directions. For the case of scalar wave scattering this doubling factor 
should be left out. The main contribution to the average transmission is given by the diffuse transport of intensity. 
This means that the two amplitudes, which make up the intensity, scatter on the same scatterers following the same 
path through the sample. In a diagrammatic language these processes are known as diffusons or ladder diagrams. An 
example of such a scattering process in drawn in Fig.|l|, in the following we will depict the diffuson by close parallel 
lines and omit the scatterers for clarity. 
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In the bulk the diffuse intensity obeys simply the diffusion equation, but the precise prefactors from the coupling 
to the outside have to be determined from the Schwarzschild-Milne equation 

C{z)^S{z)+ ( dz'M{z,z')£{z'). (4) 
Jo 

This is a self-consistent transport equation which gives the diffuse intensity C resulting from a source S. The kernel 
M describes the intensity decay between two scatterings. In the ladder approximation (without internal reflections) 
it is 



M(z,z')= / 
Jo 



The source term S is given by the single scattered incoming intensity. Our incoming plane wave leads to a source in 
the Milne equation 

_ Att 

Siz) = ntt\nf = eM-^/M). (6) 

where n is scatterer density, t is the matrix, giving the mean free path £ = 47r/(ntt). The transport equation (Q) 
now yields for an incoming diffuson 

_ iTTTliflg) L^Z 

where ti describes the limit intensity of a semi-infinite system, see Ref. 0. 

Now the prefactors are known, we continue with the diffusion equation, which holds a few mean free paths away 
from the surface. We generalize to the case where the amplitudes making up the diffuson have a non-zero momentum. 
For the corresponding diffusion propagator propagating from the point z to the point z' it holds that 

-V'C^'iz, z') - ~dlC^'\z, z') + Q^C^'iz, z') = ^^5{z - z'), (8) 

with approximate boundary conditions £'°*(0, z'\ q) = z'; q) = 0. Here Q is a two-dimensional vector, defined 

as the difference of perpendicular momenta of the two incoming amplitudes making up the diffuson. If coherent waves 
impinge with a different angle on the sample, this net transverse momentum of diffusons need not be zero and the 
diffusons decay exponentially with the inverse decay length equals \q\. The solution of the diffusion equation reads 

^int/ / „N _ 127r sinh(|Q|z<) sinh(|Q|(L ^ z>)) 
^ ' P |Q|sinh(|Q|L) ' 

with z< = min(z, z') and z> = max(z, z'). Below we also need incoming diffusons with transverse momentum, they 
are obtained by combing the prefactor of Eq. with the z> dependence of Eq. (||) . 

The total transmission is obtained by integrating over the outgoing channels. In the experimental geometry of 
Ref. [|l5|, the outgoing radiation was collected in a integrating sphere. Only outgoing diffusons where its amplitudes 
are exactly in phase (i.e. have opposite phase) are leading after this integration. Therefore the outgoing diffuson can 
have no transverse momentum. The transport equation (jj) yields an outgoing diffuson from a unit source 

^ k z , , 

^o«t ^jj^. (10) 



The total transmission is obtained by connecting the incoming to the outgoing diffuson, as was shown in Refs. [El 



{Ta) = 2 X (11) 

Here we included the doubling from the two polarization directions. Also integrating over all incoming directions 
yields the conductance 

a 

with N = 2k'^A/Aii this equals our previous definition (|l]). 
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III. CUMULANTS OF THE PROBABILITY DISTRIBUTION 



In this section we introduce the probabiUty distribution of the total transmission of scalar waves and we discuss 



some of its properties. The corrections for vector waves will be made in section VIII. We will link the moments of 



the distribution to diagrams. The moments of the probability distribution P{Ta) can be extracted as 

{T^,)^ J dTaP{Ta)Tl (13) 

In a diagrammatic approach the fc-th moment can be represented by a diagram with k diffusons on both incoming and 
outgoing side. The fc = 1 term is the average total transmission (Ta), as given by the Schwarzschild-Milne equation in 
Eg. (pi]) . This quantity is given by a single diffuson and is thus independent of channel-to-channel correlations. The 
second moment can be decomposed in the first two cumulants: 

(r2) _ {TaY + (ri)cum ^ ^ ^ ^^^2^^^ ^^^^ 



the double brackets denote cumulants normalized to the average. Diagrammatically the second moment is depicted 
in Fig. Ij. The decomposition in cumulants will prove useful as each cumulant corresponds to a different number of 
interactions between the diffusons. In the first term, Fig. ^(a) there is n o int erference; it factorizes in the average 
transmission squared (apart from a small correction discussed in section VIII ). The second term, Fig. ||(b), is the 



second cumulant ((T^)). It gives the variance of the fluctuations. Interactions between two diffusons are responsible 
for the presence of this second cumulant. The interference process is in a diagrammatical language the so-called 
Hikami-box, it is depicted as the shaded square in the figure. In the box two amplitudes of the incoming diffusons are 
interchanged, causing a correlation between the outgoing ones. Precisely the same process was studied in the more 



eneral context of the long range correlation functions 1 23 23,E3,M . In the work of De Boer et al. H and Garcia et al 



the correlation between two diffusons with different frequency was measured. This is the so called C2 correlation 
function: C2(A(jj) — {{Ta{uj)Ta{uj + Auj))). For our case find that {{T^)) — 6*2(0) and thus corresponds to the peak 
value of this correlation function. In the next section we will calculate this cumulant in detail. 

Similar to the second moment, one can distinguish three different contributions to the third moment, 

^^'^ =l + H{Tl)) + {{Tl)). (15) 



{Ta 

The corresponding leading diagrams are drawn in Fig. ^. The first term. Fig. |^(a) , again corresponds to the transmis- 
sion without interference. The second term. Fig. ^(b), is reducible in a single diffuson and a second cumulant diagram. 
From the figure it is clear that this decomposition can be done in three ways which is refiected in the prefactor of 
((T^)) in Eq. (p^). The third contribution stands for the third cumulant in the distribution and expresses the leading 
deviation from the Gaussian distribution. This is the term we are mainly interested in. It consists of two related 
diagrams: Fig. ^(c+d). The three intensities can interfere twice two by two, or the intensities can interact all three 
together, with a so called Hikami six-point vertex. Both contributions will prove to be of the same order of magnitude. 
The strength of the effect can be easily estimated using the interpretation of l/g as a interaction probability. By 
looking at the diagram, the third cumulant is proportional to the chance of two diffusons meeting twice, thus of the 
order l/g^. We find the basic result 

{{TD) « {{Tl))\ (16) 

The rest of the paper essentially consist of proving this relation and determination of the prefactor. Finally we will 
compare this relation to the experimental data of Ref . jl^] . 

Because only a limited number of channels is sampled in an experiment, the law of large numbers predicts a 
distributi on wi th some non-zero width even if we only would consider disconnected diagrams. However, we will show 
in section VIII that this effects brings only a negligible contribution to the measured cumulants. The large fluctuations 
are almost only due to the interference. 

Note that in our calculations diagrams with loops are neglected. An example of a loop diagram is the C3 (or 
universal conductance fluctuation) contribution to the second cumulant [^J|j2^ . One of the C3 diagrams is a second 
cumulant diagram where the outgoing diffusons are again input for another Hikami-box. (It is like glueing two second 
cumulant diagrams after each other.) This diagram contains two Hikami boxes, therefore it gives a contribution of 
order to the second cumulant. In general, one easily sees that in order to create a loop, one needs a higher number 
of interference vertices. Therefore, these diagrams are of higher order in 1/g, and we did not calculate them. The 
leading contributions to the cumulants are by far sufficient for the description of the experiment of Ref. . 
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IV. THE SECOND CUMULANT 



In this section we recover the results for the second cumulant. This quantity was often calculated in literature. In 
the work of the Boer et al. the frequency correlation in the total transmission was calculated using a Langevin 
approach introduced by Spivak and Zyuzin [Q, see also Pnini and Shapiro Here we obtain the same results 

using a diagrammatic technique |]6|,^,p7[ . We calculate the diagram of Fig||(b). The interaction vertex of the diffusons 
is the Hikami four point vertex or Hikami box The Hikami box is depicted in Fig.^; from the figure it is easily 
seen that the vertex interchanges two amplitudes of the incoming diffusons. As is shown in Fig.^, the full vertex is 
(in second order Born approximation) the sum of one bare vertex and two vertices with an additional scatterer. (If 
the scatterers are near resonance, the second order Born approximation is no longer valid and more diagrams are 



relevant. However, for the four point vertex it was explicitly proven in Ref. I2t 
different value of the mean free path.) The summation of the diagrams yields 



that this will only be reflected in the 



i?4 = 



Qa-qa' 



1^ 2 

qfc-qb' - ^{(ia 
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(17) 



To this vertex two incoming diffusons, a and a', and two outgoing diffusons, b and 6', are attached. The q^ denotes 
the three dimensional momentum of the corresponding diffuson. The momenta are pointing towards the vertex. For 
our second cumulant diagram of Fig.H(b) we have numbered a — 1, a' — 3, b — 2 and b' — 4. 

As stated above, only the outgoing diffusons where the amplitudes have exactly opposite phase are leading in 
the total transmission measurements. Therefore the transversal momenta of the outgoing diffusons must be zero 
'32 — Q4 — 0> but as amplitudes are exchanged in the Hikami-box, the incoming diffusons can have a transversal 
component in their momentum. Using momentum conservation one has Qi — — Q3 = Q 

After Fourier transforming Eq. (uTu in the z— coordinate we find 



Hi 



487rfc2 



Qi-Qs 



(18) 



1 = 1 



Here dz^ stands for differentiating towards the z— coordinate of the diffuson with the same number, in this case Ci{z). 
The terms inside the sum in Eq.(p^) are neglected. According to the diffuson equation, they yield approximately 
a source term near the boundary of the sample. In the integral these contributions are of the negligible order. In 
the plane wave limit of the incoming beam, all transverse momenta are absent and the diffusons are simple linear 
functions, given by Eq. (0) and (|lo|).. One obtains the known result for the second cumulant jj] 

((Tf)) = {Ta)-^ J J AxAy j\zHiC^{z)C2[z)C3{z)Ci{z) 



(19) 



Taking only this second cumulant we find for the moment a Gaussian probability distribution in the plane wave limit 



P{Ta) 



4 ^-^^ 



{Ta)f 



(20) 



Let us now study the influence of the beam profile on the correlation. If the spot of the incoming beam is 
finite, amplitudes with different transverse momenta are present. They can combine into incoming diffusons with 
perpendicular momentum. We suppose that the incoming beam has a Gaussian profile. It is decomposed into plane 
waves defined in Eq. (|^) (for convenience we assume perpendicular incidence) 



Pq 
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(21) 



where po is the beam diameter. In order to have two diffusons with a momentum Q and — Q, 
incoming amplitudes combine to a weight function 



we find that the four 



d^Pid^Pg 0(Pi)0*(Pi + Q) 0(P3)0*(P3 - Q) - e-PoQ'/\ 



(22) 
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The second cumulant is now found by first calculating the Q— dependent correlation. This is done similar to the 
calculation above, but now the Qi-Q^-teims in Eq.(|l^) should be taken into account, as well as the Q-dependence of 
the incominff diffusons, see Eq.(|). (The out going diffusons still have no such momentum.) The total cumulant is now 
found by integrating over the momentum with the corresponding weight. For a Gaussian beam profile one finds 

((^a)) = £| / e~PoQV^F,i\Q\L), (23) 

with ^2(2;) = [sinh(2a;) — 2a;] /[2x sinh^ x]. If the incoming beam is again very broad, po ^ L, only the term F2{\Q\L = 
0) = 2/3 contributes and one recovers the plane wave behavior ((T^)) ~ 2/3(7. Note that this agreement is found 
by identifying the area of a Gaussian profile with A = Trpg. This definition is somewhat arbitrary and other choices 
are also possible [ p^ . After fixing this definition, no freedom remains and we will see that for the third cumulant a 
Gaussian profile leads to other results than a plane wave. 

The second cumulant decreases as 1/ Pq at large pQ. In a real space picture it is evident that the correlation increases 
if the two incoming channels are closer to each other, i.e. if the beam diameter is smaller. In the experiment of Ref. 
I psf the focus was kept small in order to minimize the dimensionless conductance g and therefore to maximize the 
fiuct nations. 



V. THE THIRD CUMULANT 



We now discuss the calculation the third cumulant. As mentioned, there are two processes contributing. One with 
two four point vertices which we term {{T^))c and one with a six point vertex ((r^))^, where we have chosen the 
subscript according to Fig. |^. This process was, to our knowledge, not studied before. Our calculation follows the 
lines of the second cumulant calculation. 



Interference via two four point vertices 

First consider the diagram in Fig. ^(c). We have labeled to incoming diffusons with odd numbers, the outgoing ones 
with even numbers. Two incoming diffusons, Ci and £3 meet at a position z. In a Hikami box the diffusons £1 and C3 
interfere into £2 and an internal diffuson £78*- The £2 propagates out, whereas £7^* interferes again at z' with incoming 
diffuson £5 into two outgoing ones £4 and £g. Apart from this process, also three other sequences of interference are 
possible. This means that the diffusons can also be permuted as: (£1, £3, £5, £2, £4, £0) {£^,£5, £1, £4, £5, £2) — > 
(£5, £1, £3, £5, £2, £4)- We will denote the sum over these permutations as X^pcr- diagrams can also be 

complex conjugated, there is also an combinatorical factor 2 for all diagrams. (Note that complex conjugation for the 
second cumulant diagram, yields no different diagram and thus it should not be taken into account.) The expression 
for the diagram of Fig. ||(c) is now 

{{TD), = {Ta)-^2Y,A t dz dz' Hi{z)Hi{z')L^{z)C2{z)U{z)U{z')C^{z')C^{z')C}^i{z,^ (24) 



^0 

In turns out that it is useful to rewrite the form of the Hikami boxes as introduced in Eq.(^7|) into an equivalent 
expression, using momentum conservation qa + (\a' + + Qb' = 0. In real space the use of momentum conservation 
corresponds to partial integration. The Hikami box is again simplified using the fact that there are no transversal 
momentum terms, or Q terms, for the outgoing diffusons. Also the source terms of incoming and outgoing diffusons 
are neglected. Using the numbering in Fig. ^, we obtain 

Hi{z) = ^^[29., 9.. + 2a,,a,3], H,{z') - ^i!^[2a.,9,, - dl^ + qI]. (25) 

Source terms, i.e. gf— terms of the incoming and outgoing diffusons were again neglected, but the source term of the 
diffuson between the vertices is important. As is seen with the diffusion equation (|8|), it brings 

af^4"8*(z, z') + Qi4"8*(z, z') = ^5{z - z'). (26) 
The contribution from the source term (i.e. H4{z') cx ^d^^ + Qg, H4{z) (x dzidz2 + (^22^^23) is 
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7^ rL 



'Ci£2'C3i24£5£6- (27) 

Although this corresponds to a local process (just one z— coordinate is involved), it is of leading order. Together with 
the expression coming from Hi{z') proportional to dz^dz^, we find for the total contribution of the process in Fig.||(c) 



{{Ti)), = {T,r^ 



;5 .2 



) dz/:i(z)4(z)£3(2)^ dz' C'^{z')C5iz')C'6{z')dzC"'Hz,z') 
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/ dz[dz^dz., + dz^dzs + dz^dzi + dz^dz., + dz^^dz^ + dz^dz^] x 
Jo 
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Cx{z)C2{z)C^{z)U{z)C^{z)C(,{z), (28) 

where £'(z) denotes the derivative towards z of C{z). Calculated for a plane wave it gives 

22 

((Oc = Y^, (29) 
which is indeed proportional to as predicted. 

Contribution of six point vertex 

There is another diagram contributing to the third cumulant which is of the same order of magnitude as the process 
calculated above; it is depicted in Fig. ^(d). It can be thought of the following way: The use of Hikami box in the 
previous section assumes that the outgoing legs scatter at least once before they propagate out or interfere again. This 
is a reasonable assumption for the outgoing diffusons, but for the internal diffuson C™^ it also possible that coming 
from z it directly, i.e. without scattering, interferes again at z'. This process is not included in the calculation of the 
previous section but has to be studied separately. The unscattered intensity decays exponentially over one mean free 
path, therefore this process is only important if z and z' are within one mean free path. We denote the corresponding 
vertex as Hq, as six diffusons are connected to this diagram. This diagram was also already calculated by Hikami 
p8| . Also here, the dressings of the diagrams have to be added to the bare diagrams. Taking rotations of the depicted 
diagrams into account, there are 16 diagrams in second order Born approximation. It is not allowed to dress the bare 
six-point vertex (leftmost r.h.s. diagram in Fig. ^ with a scatterer which connects two opposite propagators. This 
dressing gives also a leading contribution even if the dressing performed with an arbitrary number of scatterers, but 
the resulting diagram is the same as the composed diagram with two four point vertices. Fig. |3|(c), and should thus 
not be counted. Yet this observation is useful to check the combinatorical ratio between the six point vertex and the 
composed diagram: the forbidden dressing can be performed in three ways. As the diagrams can also be complex 
conjugated, there is also a factor 2 for all diagrams. 

In the lowest order of {q€) we find for the six point vertex 

= 96^fc4 [qi-qz + q2-q3 + q3-q4 + q4-q5 + qs-qe + q6-qi + 2_.?«^]- (^o) 

i 

Hikami's original expression can be recovered from this using momentum conservation. After a Fourier-transformation 
in the z— direction the six-point vertex yields a contribution to the third cumulant 



f A 



■ / dz[dz^dz2 + dz^dzs + dzsdzi + dz^dz^ + dz^dz^ + dz^dz^] 

JQ 
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Ci{z)£2{z)C3{z)C4{z)£5{z)Ce{z). (31) 

Here we used the fact that all outgoing diffusons have zero transversal momentum. Therefore all QiQj-terms are 
absent. In the limit of an incoming plane wave we find a contribution to the third cumulant 
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{{T!)), = (32) 

The contribution from the source term, i.e. Eq.(p7|) of the previous section, exactly cancels the contribution from the 
six-point vertex. The cancellation seems plausible as one does not expect short distances properties to be important 
in the total process. Nevertheless, this cancellation depends on the precise form of the Hikami four-point vertex in 
Eq.(25). If we use other equivalent forms of the Hikami-box the contributions of the single and double integral in 
Eq.(28) are shifted with respect to each other and a full cancellation is not present. Of course, neither the result 
for Eg. ([2^) nor the final result for ((T^)) relies on this choice. The precise mechanism behind this is not clear 
to us. However, using the cancellation we only need to consider the term in Eq.(p8[), which comes from iJ4(z) cx 
(^zidz2 + ^22^231 H^lz') cx dz^dzQ and the permutations. We thus obtain for the third cumulant 

m)) = {{Ti)).+{{Ti))a 

-iTa)-'^^^Y. dz £,{z)C',{z)£,{z) dz'£',{z')£,{z')£',{z')d.C'^\z,z'). (33) 

per ^ ^ 

In the next section we calculate this expression for various incoming beam profiles. 

VI. INFLUENCE OF INCOMING BEAM PROFILE 

Now that the leading interference processes are known, inserting the diffusons gives the final value of the third 
cumulant. We first consider the simple case of incoming plane waves. As there can be no transversal momentum 
difference in the incoming amplitudes, all Qi vanish. As a result all diffusons are simple linear functions of z. And we 
find from Eq. ( |3^ ) 

{{T!)) = ^ = ^{{T^)f plane wave; po » (34) 

In practice, however, we deal with a Gaussian beam with limited spot size, influencing the cumulants in two ways. 
First, if the spot size decreases to values comparable to the sample thickness we have to convolute over a range of 
incoming momenta, just like we did when calculating the second cumulant. Second, the Gaussian profile brings an 
extra geometrical factor as we will show below. 

We need the expression when diffusons with arbitrary momentum are connected to the Hikami boxes. The outgoing 
diffusons still have no transverse momentum. Because of momentum conservation, the transversal momentum Q7 of 
the diffuson connecting the two four-boxes must equal Q5. The integration over the possible momenta results again 
in a Gaussian weight function. From the definition ( pl| ) we derive 

d^Pid^Psd^Ps 0(Pi)r (Pi + Qi) ^PsWiPs + Q3) '^(P5)0*(P5 + Q5) = e-''«(«i+'?3+Q5)/8 (35) 

Momentum conservation is used to eliminate also Q5 and reduce the integration to two transversal momenta. The 



final result for the third cumulant is obtained by inserting the momentum dependent diffusons into Eq.(33). This 
gives 



4 r 

«^a)) = Y6^ J '^''5i^''33 exp{^pl[Ql + Ql + {Q,+Q3)^]/8} F3 (|qi |L, IQ3 |L, |qi + Q3 |L) 

With 

F3(a:i,a;3,a:5) 



(36) 



(xi -I- 2:3)^X5 COSh(xi + 2:3) {xi — X3)'^X5 cosh(xi — X3) 



Yl 

{xi -I- 2:3 -I- X5)'^{xi +X3- X^y {xi - X3 + X5)'^{xi - X3- X5y 

(xi + a;3)a;5 sinh(.Ti + 2:3) (xi - a;3)a;5 sinh(xi - 2:3 



(xi -f- 


X3 -|-X5)(xi -I-X3 


- X5) 


(xi - 


X3 + X5){xi - 


- X3 - 


x^) 


{xi 4 


X3) cosh(xi -I- X3 - 


1-2x5) 




- X3) cosh(xi 


-\-X3- 


2X5) 




4(xi -f X3 + X5)2 






4(xi -1- X3 - 


x^f 




(xi - 


X3) cosh(xi — X3 - 


H2X5) 


^ (xi - 


- X3) cosh(xi 


- X3 - 


2x5) 




4(xi - X3 -1- X5)2 






4(xi - X3 - 


X5)2 





[x5 sinh(xi) sinh(x3) sinh^(x5)] , (37) 
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which is the main result in this paper. We study again the behavior if the beam diameters are wide. In the limit of 
large beam diameter (po ^ L) one finds -F'3(0, 0, 0) = y|, this means for the third cumulant ((T^)) — 4^3(0, 0, 0) /3g^, 
or 

16 

{{T!)) = y {{T^)y, Gaussian profile; po » L, (38) 

which differs by a factor | from the plane wave limit Eq.(^). This is purely a geometrical effect, depending on the 
profile of the incoming beam. In a real space picture this effect is best understood. The correlation depends on the 
distance: it is strongest if the incoming intensities are close. Therefore it is not surprising to see the influence of 
the overlap. In Rcf. this geometrical factor has been calculated for higher orders also (the area of a Gaussian 
beam is defined differently there) . For the experimental relevant case that the beam diameter is roughly equal to the 



thickness, we calculated Eqs.(p6D and (37) numerically. It then turns out that the behavior of Eq. (p8[) is actually seen 
for a large range of beam diameters. The increase of the correlation for smaller beams, turns out to be roughly the 
same for both the third cumulant and the second cumulant squared. All corrections to ( |3^ ) turn out to be relatively 
small, as we will discuss below. Apart from this advantage, errors in the sample thickness and the mean free path 
cancel by presenting the results as the ratio between the second cumulant squared and the third cumulant. 



VII. INFLUENCE OF INTERNAL REFLECTION 



In this section we calculate the influence of internal reflection on our results, and show that it is small. It was seen 
in previous work ||3l|,^ that surface reflection decreases the C2 correlation. In Eq.(^ corrections from boundary 
reflection cancel partly. We did not calculate the influence of internal reflections for the general case, but only for the 
case of very broad beams (i.e. only for Q independent diffusons). One expects that this behavior may be extrapolated 
to the Q— dependent case. At least, for the second cumulant this is a good approximation ||2^. The Q— independent 
diffusons in the presence of internal reflections are p]| 

ra ^ 4:nTi{fJ,a) L - Z + Zq ^ ^ _ k Z + Zq 

fcM/i, L + 2Z0 ' £L + 2zo' ^ ' 

here zq is the extrapolation length. In the definition of g and Tq , we replace L hy L + 2zq. If internal reflections 
are absent Zq equals 0.71£ and the corrections, which are of order zq/L, are often negligible. With internal reflection 
present zq increases and should taken into account |2ll . The correlations are known to decrease if internal reflections 
are present . In first order of zq/L the second and third cumulant behave as 

i{T!)) = ^(l-3'-^), ((r3)) = #,fl-l^^V (40) 



35 V L/ ' " 155^ V 2 L 



Therefore, the central relation (38) has a correction 

((T^)) 16 / 3 zo 



((T2))2 5 V 2 L 



(41) 



The experimental determination of the index of refraction of the sample, which determines zq, is difficult p^ . Fortu- 
nately, the correction is rather small for the experimental situation considered. 



VIII. CONTRIBUTIONS FROM DISCONNECTED DIAGRAMS 



So far the leading contributions to the second and third cumulant have been calculated. They are given by the 
connected diagrams in Fig. |^(b) and Figs. |^(c-|-d), respectively. Yet there are also contributions to the second 
and third cumulant from disconnected diagrams. The diagram in Fig. |^(a) gives an additional contribution to the 
second cumulant, and likewise, the diagrams of Figs, ^(a) and (b) give a contribution to the third cumulant. These 
disconnected diagrams correspond to cumulant contributions which are not (fully) due to interference. They describe 
effects which have little to do with the interference effects we are after. Here we calculate their contribution and show 
that they are small. 



9 



Extra contributions to second cumulant 



We first explain the contribution of Fig. |^(a) to the second cumulant. As a start we use the model of a waveguide. 
We assume that the disorder couples one incoming mode a to all outgoing modes. A waveguide has discrete modes, 
and for the moment we assume that different outgoing modes are uncorrelated. The second moment (T^) is split 
into a connected part (r^)con, Fig- §(b) and a disconnected part (T'^)dis, Fig. ||(a). The total transmission is the 
summation over all outgoing modes Ta — J2b "^ab- The disconnected part of the second moment is 

N,N N 

{Ta) dis = ^ {TabiTab^) + ^ {Tab^Tab^) 
bi^b-2 bi=b2 

= N{N~l){Tab? + N{Tl,) 

= {Ta)^+N{Tab)\ (42) 

where N is the number of modes supported by the waveguide. For the last equality sign the averaged second moment 
of the intensity speckle is given by the speckle distribution function Eq. (||), (T^j) — 2{Tab)^- From Eq.(^2|) we see that 
the disconnected diagram (T^)dis does not completely factorizc into the average squared {Ta)"^ ■ Therefore it contributes 
to the second cumulant. As shown above, the connected part of the second moment of the total transmission {Ta)con 
is proportional to L/Ni. For the sum of the disconnected and the connected contribution to the second cumulant one 
thus finds 

m)) = ((r.'))con + ((T2))di. = ^ + (43) 

which also holds for plane wave case. After this we turn to the situation of a diffusely scattering slab, with a finite 
focus of the incoming beam. The incoming beam will be broadened in the transverse direction by diffusion, changing 
the above result. To calculate the different contributions to the cumulants, for the moment the intensity distribution 
at the exit interface at transversal coordinates Ri and R2 is needed. 

The amplitudes, making up each diffuson, can propagate from to outgoing surface in different directions, respectively 
^'i, P2, P3 and P4. 

^{Ri)e'^'^' ^* {ri )e-'P'^'^{R2)e'^'^^ ** {R2)e"^^'^'- . (44) 

To get the contribution to the second moment of the total transmission, one first integrates over the transversal 
coordinates Ri and R2 to get the contribution of the whole exit interface to the intensity in a certain direction, then 
one integrates over all directions to get the total transmission. To obtain intensities the amplitudes need to be paired 
giving the following possibilities. The first possible pairing of the amplitudes is Pi = P2, P3 = Pi-, see Fig. ^(a), and 
brings 

'dPi /dP2 (*(i^i)**(fli)*(«2)«'*(fl.2)) =/'(0), (45) 



where we have defined I{q) as the transmission by a diffuson with transverse momentum q. Including the incoming 
Gaussian beam profile, it is proportional to, see Eq. (|9|) , (|2l]) 

lQle-Q'po/8 



Integrating over all outgoing directions results in 



(47) 



This is the just the factorizing contribution, and hence does not contribute to the second cumulant. 

The second possible pairing of the amplitudes. Pi = P4, P2 — P3, does give a contribution to the second cumulant. 
It is the diagram in Fig. ^(b). 

J dRidR2 (*(i?i)**(Pi)*(i?.2)**(fl2))G*^i(^i-^^)e-''^^(^i-^^) = /(Pi - P2)/(P2 - Pi). (48) 
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Subsequent integration over the outgoing directions Pi and P2 yields 

{TD^^el d'Ql^Q). (49) 

J\Q\<k 

The integral can be extended to infinity since I{q) is an exponentially decaying function. The contribution of the 
disconnected diagram Fig. ^(b) to the second cumulant is thus 

iTk'r d'aPia) 1 

In analogy to the waveguide, this result describes irreducible contributions from disconnected diagrams. It can be 
interpreted as the inverse of the number of independent speckle spots in transmission at the exit interface [p3^,f5| . 



Extra contributions to the third cumulant 



We apply the same method for contributions to the third cumulant. Following the waveguide argument as above, 
we find the contributions to the third cumulant of the diagrams of Figs. ||(c+d), ||(b), and ||(a), respectively, 

. . . 3 L'^ Q L 2 , 
((^a))= 5^+2]^ + ]^- (51) 

The first r.h.s. term is the connected diagram. It is clear that the second term, the diagram Fig. ^(b), gives a much 
larger contribution to the third cumulant than the third term, the diagram of Fig. §(a), as the diagram Fig. ^b) is 
already enhanced by some interference. In the following we only consider this diagram. 

As can be seen from Fig. ^(b) there are three possibilities to combine the three diffusons into two connected diffusons. 
and a single diffuson. Attaching outgoing directions to the amplitudes at the exit interface gives (see Fig. (^, 

'^iRi)e'^'^'^^Ri)e-'^'''^'^{R2)e'^''^''^*{R2)e-'^'"'-^iR3)e^^ (52) 

It is clear that there are six possibilities to pair the outgoing directions into intensities. Integrating over the transversal 
coordinates Ri, R2 and ^3 gives the following six contributions, 

/con(0)/co„(0)/(0) /con(Pl " P3)/con(P3 " Pl)I{0) 

Icon{0)Icon{Pl " P5)I{P5 " ^^l) ^con (0)/con (^3 " P5)IiP5 ~ Pz) 

IconiPl - P5)IconiP3 " Pl)IiP5 " ^=3) /con(^'l - P3)IconiP3 " P5)I{P5 " ^^l), (53) 

where Icon denotes the transmitted intensity coming from the connected part of diagram Fig. ^(b). The third cumulant 
is in the discrete mode model given by 



{{t!)) 



(Ta) 



b,b'.b" b b,b' 



(54) 



Similarly to the second cumulant we insert all pairings of Eq. ( p^ ) and find for the contribution from disconnected 
diagrams 

((r„3))dis = J7FT3 I <i^Qld^Q2d^Q3 [/con(0)/con(Ql)/(-Ql) + /co„(Ql)/con(Q2)^(Q3)] • (55) 



(Ta) 

The integrand is dominated by its first term, which is depicted in Fig. |^. We explicitly evaluated this term. The 
calculation follows completely the line of the second cumulant calculation, yet it is slightly more complicated as it has 
a Q— dependent outgoing diffuson. The integration over Q2 and Q3 gives a factor (tt/c^)^. Comparing the result with 
the last two terms in third cumulant in the waveguide model, Eq. (pl|), we define N* as 

((T,^)>dis = 6((r2))_/Ar*. (56) 

The number N* is inversely proportional to the contribution of disconnected diagrams to the third cumulant for a 
Gaussian profile. 
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Polarization effects 



The vector character of the hght has not been taken into account yet. The two independent polarizations of each 
outgoing direction effectively double the number of independent speckle spots N. Thus for an incoming plane wave 
of unit intensity with fixed polarization the total transmission and the conductance g are twice as large as they 
would be in the scalar case. As we work with normalized cumulants, this effect reduces only the value of the second 
cumulant (oc 1/g) and the value of the third cumulant (cx l/g^). Therefore, it is immediately seen that the relation 
((T^)) oc ((T^))^ (Eq. ^) is not affected. The vector character does reduce the correction Eq. ( |50| ) by a factor 2. For 
the experimental data of Rcf. [|l5| the number of modes N as well as N* is listed in Table | (including the doubling). 

Summarizing the previous sections we have included three corrections. We first obtained the result for very broad 
Gaussian beams, in the large L/£ limit, Eq. (^8|). The first correction was the influence of a finite beam diameter, it 
changes the diffuse intensity from linear into an exponentially decaying, see Eq. (^. This correction is contained in 
Eq. ( p7|) and Eq. (|36|). The presence of internal refiections also changes the spatial dependence of the diffuse intensity 
resulting in a correction Eq. (p]). The third correction is of another nature, it is the only process which does not 
come from interference, but from disconnected diagrams. Only this term depends the number of modes, which in the 
vector case is twice as large as in the scalar case. 

IX. COMPARISON WITH EXPERIMENTS 

The data set found experimentally in Ref. is reproduced in Table | and Fig. ||. The experiments reported there 
were performed with seven different samples. The experimental setup and measurement technique used is extensively 
described in Ref. Samples consisted of 36 vol.% rutile Ti02 pigment on a transparent substrate. The extrapolation 
length was estimated from the effective index of refraction to be zq « 1.1/im. The absorption length £a was determined 
to be ~ 70/im. Different values of the conductance g were probed by taking various sample thicknesses and by varying 
the beam diameter. The fluctuations in the total transmission were measured by varying the wavelength of the light. 

For a very broad beam we found the simple relation (^) between second and third cumulant. A weighted least 
square fit to 

{{T!)) = const.{{T^)r (57) 

of the raw experimental data yields a prefactor 2.9 ± 0.6. However, as discussed above, there are three corrections to 
be made. First, if the beam width becomes comparable to sample thickness the integrals (|3|) and ( |36| ) have to be 
performed: If the beam width reduces, g decreases accordingly and both cumulants increase in absolute size. Yet the 
precise increase is somewhat different resulting in a somewhat smaller prefactor in Eq. (|38|). We corrected each data 
point individually for its finite focus, mapping it to the infinite focus case. The third cumulant was multiplied by a 
factor which ranged from 1.03 to 1.13, as L/po ranged from 0.41 to 6.3, see Table |[ This is the largest correction, 
it changes the prefactor some 10%. Secondly, we corrected for internal reflections according to Eq.([ll|). The third 
correction comes from the disconnect diagrams. The contributions from the disconnected diagrams are substracted 
from the measured cumulants. After all these corrections the data should again obey the law: {{Tj^)) = 3.2 ((T^))^. 
The results are plotted in Fig. |^, where the points are the corrected data points and the line is the theoretical 
prediction. A least square fit gives 

((T3))^(3.3±0.6)((T2))2. (58) 

Note that there is no adjustable parameter. We find that there is good agreement between experiment and theory. 
All corrections are minor as compared to the error in the data, fitting the raw data is also in agreement with the 
theoretical value of 3.2. We recall that the major shift between fits of raw and corrected data comes from the beam 
diameter coreection. Inspecting the figure one might be tempted to make a linear fit, but in Ref. jist it was shown 
that this fit is statistically very improbable. 

X. DISCUSSION 

We have calculated the second and third cumulant of the distribution of the total transmission underlying the 
conclusions of Ref. jist and compared it to the experimental data. Both cumulants are a consequence of interference 
between diffuse channels. They were calculated with a diagrammatic technique. The inverse dimensionless conduc- 
tance, interpreted as an interference probability, is a perturbation parameter in the theory. The third cumulant is 
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proportional to the second cumulant squared. We also found a non-trivial dependence on the profile of the incoming 
beam used. The cumulants were calculated for arbitrary beam diameter, but the influence of a finite focus on the ratio 
is rather weak. Also boundary reflections were included. Our calculations confirm that the main contributions come 
from diagrams with interference processes, i.e. connected diagrams, as we have shown that the contributions from 
disconnected diagrams is small. The experimentally found ratio of the third cumulant versus the second cumulant 
squared is well described by our theory. 

The extension of the calculations to higher cumulants is straightforward. The n-th cumulant will contain (n — 1) 
Hikanii four-point vertices. So the contribution is estimated to be: {{T^}) oc Also corrections and cancellations 

from higher order vertices are present, but it is clear that the calculation becomes very laborious at large n. Recently 
two of the authors discovered that all the cumulants of the distribution function can be mapped onto the moments 
of the eigenvalue distribution of the transmission matrix p^ . The eigenvalue distribution is bimodal and was first 
calculated using random matrix techniques ]34| ] but recently its validity beyond quasi-lD was proven |35|. As the 
eigenvalue distribution is known, the entire distribution of the total transmission was calculated in the limit of broad 
beams. These results agree with calculations presented here for the first three cumulants. The experimental data 
thus also proof the first few moments of the eigenvalue distribution function. As only three moments are known, 
it is impossible to reconstruct the full eigenvalue distribution from the experimental data. The present calculation 
leads us to assume that the eigenvalue distribution in a diagrammatic approach is also given by loopless connected 
diagrams. As the ratio of first few cumulants does not depend sensitively on the beam diameter, the results of Ref. 
|16| are probably also valid in the regime where the beam diameter becomes comparable to the sample thickness. 
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TABLE I. Sample thickness, beam width, second cumulant, third cumulant for the different samples as taken from Ref. [ [l5| . 
Next are the beam diameter correction factor on the third cumulant (see section ^), and the number of modes A'^ and A'^*. 
Together they give the corrected (plotted) cumulants (last two columns). 



Sample 


Beam 


Second 


Third 


Beam 


Number 


Number 


Corrected 


Corrected 


thickness L 


diameter po 


cumulant 


cumulant 


diameter 


of modes 


of modes 


2e cum. 


3th cum. 


(in nm) 


(in ^m) 


(xlO"") 


(xlO"^) 


correction 


N 


N* 


(xlO"*) 


(xlO"'') 


30 


77 


0.36±0.01 


0.014 ± 0.035 


1.03 


388000 


300000 


0.33 


0.007 


12 


26 


0.97 ±0.03 


-0.03 ±0.25 


1.04 


46800 


36800 


0.75 


-0.18 


22 


32 


1.24 ±0.04 


0.68 ±0.28 


1.06 


88600 


72100 


1.13 


0.58 


30 


33 


1.57 ±0.04 


1.30 ±0.46 


1.07 


119000 


99300 


1.49 


1.21 


53 


35 


1.80 ±0.03 


0.91 ±0.53 


1.10 


241000 


212000 


1.76 


0.86 


30 


26 


1.90 ±0.03 


0.92 ±0.56 


1.09 


94900 


81300 


1.79 


0.78 


45 


33 


1.90 ±0.05 


1.33 ±0.43 


1.10 


187000 


163000 


1.84 


1.26 


53 


26 


2.18 ±0.03 


1.77 ±0.59 


1.10 


208000 


189000 


2.13 


1.70 


170 


27 


2.69 ±0.06 


2.02 ±0.82 


1.13 


1420000 


1430000 


2.68 


2.00 


78 


28 


2.74 ±0.03 


2.43 ±0.62 


1.11 


396000 


372000 


2.71 


2.39 


30 


17 


4.82 ±0.10 


9.1 ±3.3 


1.11 


71900 


64400 


4.68 


8.60 


30 


10 


8.01 ±0.36 


5.3 ±6.4 


1.11 


60200 


56400 


7.84 


4.47 
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FIG. 1. Left: an example of an actual scattering process; a retarded (full line) and an advanced amplitude (dashed line) 
come from the left and share the same path through the sample. Right: schematic representation of the average process, the 
diffuson (for clarity the scatterings are not drawn). 

FIG. 2. The two contributions to the second moment of the total transmission. In diagram (a) the transmission channels 
are independent; this process is of order unity and is almost completely reducible to the mean value squared. Diagram (b) 
corresponds to two interfering channels; this is the second cumulant, it is of order l/g. The close parallel lines are diffusons; 
the shaded square denotes the Hikami four point vertex. 

FIG. 3. The three contributions to the third moment of the total transmission. Diagram a) corresponds to independent 
transmission channels; it is of order 1. In diagram b) there is correlation but it can be decomposed into the second cumulant; 
this is of order . The diagrams c) and d) are the contributions to the third cumulant, 0{g~^). 

FIG. 4. The Hikami four point vertex. It describes the exchange of amplitudes of two incoming diffusons a and a' into two 
outgoing diffusons b and b' . The dots linked with the dashed line denote the dressing with an extra scatterer. 

FIG. 5. Diagrams contributing to the interaction of six diffusons: H^. To 1,3,5 the incoming diffusons are connected, to 
2,4,6 the outgoing ones. Possible rotations of the three rightmost diagrams are not drawn, in total there are sixteen diagrams. 

FIG. 6. The disconnected contribution to the second moment. In (a) the two transmissions factorize into the average value 
squared, and does not contribute to the second cumulant. Diagram (b) is much smaller, but gives a contribution to the second 
cumulant. The amplitudes making up the two diffusons propagate out in different directions. 



FIG. 7. The leading disconnected contribution to the third cumulant, the box symbolises again the Hikami vertex. 

FIG. 8. The third cumulant plotted against the second cumulant. The points are the corrected experimental data (see text). 
The line is the theoretical prediction, {(Tq)) = 3.2((T„)); no free parameters were introduced. 
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